Effective interactions between inclusions in complex fluids driven out of equilibrium 
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The concept of fluctuation-induced effective interactions is extended to systems driven out of 
equilibrium. We compute the forces experienced by macroscopic objects immersed in a soft material 
driven by external shaking sources. We show that, in contrast with equilibrium Casimir forces 
induced by thermal fluctuations, their sign, range and amplitude depends on specifics of the shaking 
and can thus be tuned. We also comment upon the dispersion of these shaking-induced forces, and 
discuss their potential application to phase ordering in soft-materials. 
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I. INTRODUCTION 

A prominent issue in soft condensed matter physics 
is the understanding of the equilibrium phase behavior 
of mesoscopic particles immersed in complex fluids: col- 
loidal suspensions Q , liquid droplets in liquid crystals , 
inclusions in lipid membranes 0], charged particles in 
electrolytes 0,... A standard and fruitful procedure to 
describe the large scale properties of inclusions in a soft 
medium consists in integrating out the numerous "sol- 
vent" degrees of freedom via ensemble averaging 0, . 
The interaction between the embedded particles is then 
described by "effective potentials" . The latter can mod- 
ify the genuine interparticlcs interactions (Coulomb, van 
der Waals,...) or give rise to entirely new effects. For in- 
stance, in some (equilibrium) cases, the external objects 
do not modify the ground state energy of the medium 
but only alter its thermal-fluctuations spectrum. The re- 
sulting entropic effective interaction has consequently an 
amplitude proportional to the thermal energy and a 
range comparable to that of the correlations of the fluc- 
tuations of the medium Q ■ Such fluctuation- induced 
(i.e. entropic) interactions are commonly referred to as 
thermal Casimir interactions in analogy with their fa- 
mous quantum equivalent 

More recently, many experimental studies have re- 
ported the organization of particles embedded in fluids 
when the latter are driven out of e quili brium by the appli- 
cation of external fields: [9lll0llTllll2l|. Extensions of the 
concepts of "effective potential" and "entropic forces" to 
out of equilibrium situations have been scarce [T^IT^Il^ . 
Indeed effective potentials can not be simply derived from 
a free energy in an out-of-equilibrium context, and only 
the instantaneous force acting on the external objects 
for a given configuration of the medium can be defined. 
The effective interactions between the host objects should 
then strongly depend on the dynamics ruling the tempo- 
ral evolution of the medium (in contrast to the equilib- 
rium Casimir case). 

In this paper, we attempt to extend the paradigm of 
Casimir effective interactions to objects immersed in soft 
systems "shaken" by external energy sources that can 



FIG. 1: Sketch of the three shaking processes discussed in the 
paper. (A): A point-like monochromatic shaking source. (B): 
Uniform monochromatic shaking. (C): Spatially uncorrelated 
colored noise. 



not be a priori modeled by heat baths. We use a model, 
detailed in section II, where both the medium and the 
objects are very simple, (i) The hosting medium is de- 
scribed by a scalar free field (j) living in a d dimensional 
space (d = 1,2,3): for instance, (j) could represent the 
height profile of a fluid interface {d = 2) or a contact 
line (d= 1) El, or the angular deviation of the director 
of a nematic liquid crystal {d = 3) ^3- (ii) The ex- 
ternal objects are two identical rigid parallel plates that 
enforce a zero field on their surfaces (Dirichlet bound- 
ary conditions). Among the numerous possible choices 
for the shaking sources, we explore three different cases 
of experimental relevance (c.f. Fig. QJ: (A) a localized 
monochromatic shaking, (B) an uniform monochromatic 
shaking, (C) a spatially uncorrelated colored noise. All 
these shaking induce effective interactions between the 
plates, which we characterize in section III by the cor- 
responding average forces, discussing their sign, ampli- 
tude and range. We compare these features with those 
of the usual thermal Casimir effect. The detail of the cal- 
culations and a more thorough analysis of the shaking- 
induced interactions are gathered in section IV. This 
more technical section also includes some comments on 
the time dependence and fluctuations of these shaking- 
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induced forces. Section V ends the paper with a synthetic 
summary of our main results and an outlook on possible 
applications. 



II. MODEL AND NOTATIONS. 
A. Energetics and Force on the plates : 

We consider two plates in a soft medium separated by 
a distance L much shorter than their lateral extension 
L|| as depicted in Fig. |2 Their thickness plays no role 
in all that follows and will be set to 0. The soft medium 
is modeled by a scalar field (j) associated to the elastic 
Hamiltonian : 



^^2 



V<^(i?) 



(1) 



with K, the elastic modulus and R = (x, r) where x is the 
coordinate normal to the plate, see Fig. El is taken 
dimensionless. 

We restrict ourselves to strong plates-field interactions, 
modeled by Dirichlet boundary conditions (DBC) on the 
plates: 0(a; = 0, r) = 0(x = L, r) = 0. 

We first focus our attention on the force that the 
soft medium exerts on the right hand side of the plate 
located at a; = 0. For a given configuration of the elastic 
field, is given by the integral over r of the stress ten- 
sor component Txx{'^^) jlS- The components Tij of this 
stress tensor can be obtained by the virtual displacement 
method taking into account the DBC 0. In the present 
geometry, Txx{r) = -^k [5j;(/)(0+, r)] , so that 



dv [a,0(O+,r)]2 



(2) 



Note that this force pushes the plate toward the negative 
X direction whatever the field conformation. 

The net force F on the plate is the algebraic sum of 
the contributions of the two sides: F = F> -|- F<. The 
force F< is given by a formula similar to Eq. (j^J, so that: 



- /dr [9,0(0-, r)]2. (3) 
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FIG. 2: Two parallel plates, perpendicular to the x axis, sep- 
arated by a distance L shorter than their lateral extension 
Ln. 



Here 7 is a generalized friction coefficient. 

A description of the dynamics of the plates is beyond 
the scope of this paper, and we assume that they are 
fixed. L is thus a constant. After partial Fourier trans- 
form with respect to t and r, equations (|4I5|) can be recast 
into: 



Jo 



(6) 



where the response function TZ^ corresponds to the dif- 
fusion kernel in the slab geometry with DBC and with 
a diffusion constant k/^. The Fourier transforms used 
in this paper are defined by /q = / dr /(r) exp(iq.x) and 
= Jdtg{t)exp{iujt). 

A generalization of our results obtained with H4I5|) 
to other slow dynamics that obey an equation of the 
form JB)) with a different kernel will be commented upon 
throughout the text. This generalization is simple pro- 
vided the dynamic linear response of (j) relates spatial and 
temporal scales through an algebraic relation l^^ ~ uj~^^^ , 
where is the spatial extension of the elastic distortion 
resulting from a localized periodic shaking of pulsation 
Lu, and z is the so-called dynamical exponent. The dif- 
fusive model described above by Eq. H4I5|I indeed fits in 
this picture with z = 2 and: 



L 



-1/2 



(7) 



B. Dynamics : 

For sake of clarity all calculations are performed within 
the simplest local and strongly dissipative dynamical 
scheme: the (j) field evolves under the application of an 
external "shaking" source rj as : 



79t0 = — TT +V 



</.(0,r,i) 



5n 

(j>{L,r,t) = 



(4) 
(5) 



III. AVERAGE FORCES INDUCED BY THREE 
KINDS OF SHAKING 

To set a reference for further comparison, we first re- 
call the expression of the average thermal force -Fcasimir 
on the plate at a; = if the whole medium is thermally 
excited by a heat bath imposing a temperature T |22l |. 
The simplest derivation consists in computing the total 
free energy of the system for a given interplate distance 
L, and taking its derivative with respect to L. An al- 
ternative approach is to consider the model Langevin 
dynamics for the field given by Ea. (|4l5|) . with 77 a ther- 
mal Gaussian white noise of zero mean and of variance 
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{T]{R,t)T]{R\t')) =2jkBTS{R- R')5{t-t'). Averaging 
over realizations of the noise leads to the same universal 
expression for the net thermal Casimir force 0: 

J^Casimir = k^T , (8) 

with Ad ^ (d- l)r{d/2)C{d)/{4:7rY/^. This universal 
mean force describes a long range attraction that depends 
on the temperature of the heat bath but not on specific 
material properties such as k and 7. 

We now calculate the average forces generated by the 
three non-thermal shaking processes mentioned in the 
introduction, with a field (jj that evolves according to 
Eq. 

A. Average force induced by a localized 
monochromatic shaking 

We start with the simplest choice for the shaking in 
Eq. Q , namely a point- like monochromatic shaker oscil- 
lating at a pulsation cug, located at an equal distance L/2 
from the two plates: 

r]{R, t) = rjsa'^ cos{LJst)S{x - L/2)S{r), (9) 

The noise amplitude has been arbitrarily split into the 
product of a microscopic volume a*^ where the medium 
is shaken, times a shaking amplitude % that has the di- 
mension of an energy density. 

Such a process could be, for example, the result of the 
forcing of a microscopic " active" particle located between 
the two plates (e.g. a magnetic particle under the influ- 
ence of an external oscillating magnetic field), in which 
case a is typically fixed whereas 77s and tOs can be exter- 
nally tuned. 

The net force F felt by the plate at a; = 0, reduces 
to its right hand side contribution since the outer 
medium remains undistorted < 0) = 0). Moreover 
the geometry is here totally symmetrical so that the force 
felt by the plate at x = L is: —F. Denoting F^ the time 
averaged force for this process, and anticipating on the 
detailed discussion of section IV-A, we use Eq. (|23|l to 
write: 

Again l^j^ = [{j / k)uJs]~ ^ is the dynamical length as- 
sociated with the shaking pulsation ujs. The adimen- 
sional function Ta(1/2,L/^„J is plotted in Fig. and 
its general expression (including the shaker position de- 
pendence) is given below in Eq. (|23|l . 

Two different limit regimes can be identified: If 
L/lcjs <C 1, the shaking period is much larger than the re- 
laxation of any interplate field excitation, so that the elas- 
tic deformations generated at x = L/2 propagate along 
the X axis up to the plates surfaces. In this quasistatic 




FIG. 3: Localized monochromatic shaking. Log-linear plot 
of the scaling function T a{1/2, L/l^^) against L/l^i^ = 
/k^^'^ujI^^ L, for various space dimensions: d — S full line, 
d = 2 dashed line, and d = 1 dotted line. 

limit, the time averaged force F^ felt by the plate at 
a; = is comparable to the elastic force on the plates un- 
der the application of a constant static "effective" pertur- 
bation ■^T]sa'^6{x — L/2)(5(r). The latter induces static 

deformations of lateral extension L on the plate. Thus, 
F^ scales as r/^a^'^ / {kL'^~^) independently of the shaking 
period. In contrast, the dynamic length is much smaller 
than the plate-shaker distance if L/lui^ ^ 1. The elastic 
deformations are exponentially screened before reaching 
the plates surfaces. Therefore, the force amplitude de- 
cays exponentially as shown in Fig. |3| 

In summary, a localized periodic shaking induces an ef- 
fective interplate interaction. Whereas J^casimir is univer- 
sal and leads to an attraction between the two plates, the 
interaction obtained here is repulsive and docs depend 
both on the elastic modulus of the surrounding "fluid" 
and on its dynamics (through 1^:^). Finally, both the 
range and the amplitude of the present shaking-induced 
force can be externally tuned by varying the pulsation lUs 
and the amplitude ?7s of the shaking. The dependence of 
the force on the shaker position as well as its fluctuations 
will be discussed in section IV. 



B. Average force induced by a uniform 
monochromatic shaking 

We now suppose that the elastic medium is everywhere 
periodically and homogeneously shaken: 

■q{R,t) ^ 7]sC0s{uJst)- (11) 

Posible experimental realizations include: the homoge- 
neous and periodic shaking of a metallic elastic line us- 
ing a magnetic field |17| | or the driving of an interface 
between two dielectic liquids using AC electric fields . 
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FIG. 4: Uniform monochromatic shaking. Linear plot of FIG. 5: Colored noise. Log-log plot of the ratio Tc(-£///f2s,)/Ad 
the scaling function Tb(L/LJ as a function of L/L, = versus L/ln, = {-y / k^^^^V^ L. The full line corresponds to 
{y/Ky^^ujl^'^L . d = 3, the dashed line to d = 2 and the dotted line to d = 1. 



One may also use flexoelectric effects in liquid crystal to 
generate such a shaking process . 

By symmetry the two plates feel opposite instanta- 
neous forces. The exact expression for the time averaged 
force Fb on the a:: = plate is computed in section IV: 



2 rrf-l 



-Fb 



■ It Tr ( 



(12) 



The dimensionless quantity Tb, see Eq. (|30|) . is plotted in 
Figure^ it is non-monotonic and changes sign periodi- 
cally. Given the plates geometry and the shaking transla- 
tional invariance, the net normal stress FB/ij|^^ depends 
neither on the plates lateral extension L|| nor on the di- 
mension d. The fast and slow shaking asymptotics can 
thus be qualitatively inferred from dimensional analysis. 
If the medium is slowly shaken {L/l^ ^ 1), the net force 
amplitude diverges as ~ 1/t^s- First, notice that a 
constant uniform source 77s between the plates induces an 
homogenous normal stress on the inner side of the plate 
(x = 0+) that scales necessarily as ~ ?7s k^^L^. 

Second, on the left hand side of the plate (x < 0), l^^ 
is the only length scale available to construct the normal 
stress /Ly~"^ ~ 'nl'^~^^t^- I^i tl^is slow shaking regime 
this latter contribution dominates, so that Tb plateaus 
to a finite value. Turning now to the L/l^ ^ 1 limit, 
if the interplate distance goes to infinity, the two sides 
of the plate face identical semi infinite media. Conse- 
quently, they arc pushed in opposite directions with the 
same amplitude and thus feel on average no net force. 
So, in this limit L/l^ 3> 1 that also corresponds to fast 
shaking at fixed L, Fb decays to zero, as observed in Fig. 



a. 

Returning to the most remarkable features, Fig. 2]and 
Eq. 112|) show that the fluctuations of a uniformally 
shaken medium induce effective interplate forces that 
strongly differ from -Fcasimir • They arc indeed completely 



tunable unlike their equilibrium analog. Their sign, am- 
plitude and range can all be controlled externally by vary- 
ing the shaking period and amplitude. 



Average force induced by a stochastic shaking: 
colored noise 



Finally the case of stochastic shaking sources is con- 
sidered. As mentioned above, if 77 is the white Langevin 
noise modeling the coupling to a heat bath, the result- 
ing field fiuctuations induce thermal Casimir interactions 
between the plates. Any other noisy shaking drives the 
system out of equilibrium. Numerous examples of elastic 
materials shaken by "active" stochastic processes are pro- 
vided by biological systems, e.g. active membranes [ll|, 
actin-myosin gels |21|, etc. 

Without attempting to describe accurately the 
specifics of the noise in a given system, we focus here on 
the simple case of a spatially uncorrelated colored noise 
of zero mean, seeking for the infiuence of its color on the 
interplate force generated. We describe the fluctuations 
of this noise by: 

(77(i?, t)iiR', t')) = rila-'Sit - t')5{R - i?'), (13) 

where the correlation S{t) is deflned by its Fourier trans- 
form also referred to as the noise power spectrum: 



(14) 



In Eq. (|13|l a is a microscopic distance chosen so that 
the mean squared noise satisfies {rj^{R,t)) = r^. Again 
77s has the dimension of an energy density. The white 
noise limit corresponds to fig cx), with rila'^Vl~^ kept 
constant. 
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The mean net force Fc on the x 
computed in section IV-C: 



plate is exactly 



2 d Lf-^ 



(15) 



The dimensionless ratio Tc/A^ is plotted in Fig. [SJ us- 
ing Eq. (123 ■ We recall that Ad measures the am- 
plitude of the thermal Casimir force {Ad = [d — 

i)r(d/2)c(d)/(4^)'^/2). 

Postponing a more involved discussion to next section, 
we emphasize here three points. First, the color of the 
noise has no influence on the sign of the force: the two 
plates attract each other whatever llg and 77s. Second, 
Tc decreases monotonically with L /l^^ , in other words 
the larger the noise correlation time fl^^, the stronger 
the force felt by the two plates. Third, Fig. jsj indicates 
that Tc{u)/Ad decays as l/(2u^) when u ^ 1. This 
implies that the attraction is long-ranged and decays only 
algebraically with the interplatc distance. A connection 
with the Casimir result appears in this regime: if one 
sets rj^a'^fl'^ = then in the limit white noise 

limit Q.S Qo, Fq converges as expected to -Fcasimir = 
AdfcBTLff-VL'*. 

To conclude, we have shown that Casimir-like attrac- 
tions can be induced by colored noisy shaking. These 
attractions decay algebraically at long distance as in the 
thermal case and contrary to the two former monochro- 
matic cases. Their amplitude depend explicitly on the 
material properties and their scaling on the noise color. 



IV. SHAKING-INDUCED FORCES : DETAILED 
ANALYSIS. 

In this section we provide a more detailed analysis of 
the interactions induced by the three kind of shakings: in 
addition to the explicit derivation of the average of the 
corresponding forces, we also comment on their temporal 
fluctuations. 

We start with the derivation of the tool that will al- 
low us to compute the forces on the x ~ Q plate in 
all situations, namely the linear response function dxTi- 
of the field's gradient at the surface of the plate to a 
point-like shaker in the medium (by definition, the ker- 
nel TZ is the solution of the dynamic equations (|4I5|I . with 
i^{R,t) = 5{R - R')5{t - t')). We compute this quan- 
tity using an image method which conveys an intuitive 
picture |23 |: as in electrostatics, these equations can be 
solved replacing the constraint at the boundaries lO by 
the addition of an ad hoc distribution of image sources 
outside the integration domain. 

We consider first a shaking source located in the in- 
terplate region Q < x' < L and compute the response 
dxTi-^ on the a; = 0+ side of the plate. In the present 
slab geometry, the images are the mirror images of the 
original source through virtual reflecting planes located 
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FIG. 6: Distribution of the image sources for the two plates 
geometry and Dirichlet boundary conditions. The filled circle 
is the original source and the dashed circles its mirror images 



at X = nL as depicted in Fig. |S1 Consequently, the re- 
sponse 7?,^q(x,x') can be expressed as an infinite sum 
over all the images contributions: 



T^w.q {2nL + X — x) 
7?-t^,q {2nL — x' — x) ■ 



(16) 



where TZ the diffusion kernel in an infinite domain 

exp {^-yjq^ + i \x — x'\j 



'^LJ,q('^i ) 



2k-\/ q"^ +ilu 



(17) 



with Vl — 1. Performing the sum H16|l and taking the 
derivative with respect to x yields: 



sinh 


^iqL)^+iiL/L)^l~x 


/L)_ 


Ksinh 







9,7e>q(0,x') 



(18) 

A simpler calculation can be carried out if the source 
is in the negative x region {x' < 0). In this case the only 
image source is the symmetric in the a; = plane of the 
original source, so that the response of the field gradient 
on the X = 0~ side of the plate is described by: 



a,7^< (o,x') = 



(19) 



A. Localized Periodic Shaking. 

If the only driving is a single monochromatic shaker 
located between the plates 3± x = x^^ 



r] ^ iisa'^S{x - Xs)d{r) coa{LUst), 



(20) 



then the general expression of the net force acting on the 
X = plate is simply obtained from the definition of the 
response function dxTZ^'- 



d\2 



|c>,7^>,q(0,X,)| 



[1 + C0S{2 LOst)] 



(21) 
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sus Xs/L. This corresponds to case A in the quasi-static limit 
L/lus "C 1. The full line corresponds to d = 3, the dashed 
one to d = 2 and the dotted line to d = 1 
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FIG. 8: Log-linear plot of the scaling function Ta- The full 
line corresponds to d = 3, the dashed line to d = 2 and the 
dotted line to d = 1 



with the short hand notation J^ = J dq/{2ny~^ [2J|. 

There is obviously no force on the left side of the plate, 
so that using H18|l we obtain the exact expression for the 
instantaneous force on this plate. Its time average 
describes a repulsion between the plates: 



d\2 



Xs L 



(22) 



with, 



Ta(u,w) - 



sinh 



\/ + iv"^ (1 — It) 



sinh( \J + iv'^) 



(23) 



We now give a closer look at the asymptotics corre- 
sponding to the fast and slow shaking limits. 

Quasi-static shaking: We first consider the regime 
where the driving period is much larger than the relax- 
ation time of any interplatc field deformation, i.e. the 
dynamical length is much larger than the intcrplate dis- 
tance {L/l^^ 1). 

The limiting value TA(it, 0) can be exactly com- 
puted by integrating over the in plane q modes 
before summing up all the contributions from 
all images in Eqs. l)16li>l|) . If = 3 Ta(w,0) = 
du [tt (1 - u) cot (ttm) -2Tu + 2{l-u)ip (u)] /(327r), with 
^ the Euler ■i/'-function, and T the Euler constant [25l |. 
if d = 2, Ta(m,0) = [l + 7r(l-u) cot(7ru)] /(Stt), and 
if d = 1, Ta('u, 0) = (1 — u)^/4. The scaling function 
Ta{xs/L,0) is plotted in Figure [7| (log-log plot). The 
fast decay for values of Xg/L larger than i witnesses 
that the elastic distortions are essentially occurring in 
the L/2 < X < L region. Conversely a shaker close to 
X = induces strong variations of the elastic field in the 
vicinity of the plate, leading to an algebraic divergence 



of the force when x^/L (if d > 1). 
The force dependence on the shaker position, can also 
be qualitatively understood thanks to an electrostatic 
analogy. In the present quasi-static limit, the kernel TZ^ 
corresponds to the inverse of the Laplace operator (with 
DBC on the plates). Consequently cj) plays the role of 
an electrostatic potential, the plates mimic a grounded 
planar capacitance, and the shaking source a point like 
particle bearing a charge rjsa'^. The force on the plate is 
then analogous to the electrostatic force produced by a 
point like charge |26l |. From the image expansion H16() 
we know that this force is identical to the one produced 
on a virtual plane at x = by an infinite number of 
(+77sa'^) charges located at a; = 2nL + Xs and of (— TysO"^) 
charges at a: = 2nL — Xg- Hence: 

-If Xs < L (c.f Fig. El): the electric field on the plate 
surface is dominated by the image charges at = Xs, so 
F ^ —{r)sa'^)^/{KXg~^). Higher order reflections produce 
image sources that form dipoles viewed from x = 0. 
Their subdominant contributions to the electrostatic 
force actually reduce the above estimation. This implies 
that the further the second plate, the bigger the force 
produced by a single shaker. 

-if Xs ^ L (c.f Fig. 10): the image charges form dipoles 
located at x = (2n -I- 1)L. The electric field on the plate 
is dominated by the influence of the dipole located at 
I a; I ~ L (the closest from the a; = plate). The force 
therefore scales F^ {Vsa'^)^{l - Xs) Lf / [kL'^-^) . 



Fast shaking: If Ljl^^ <C 1, the electrostatic analogy 
only provides a partial description of the force generation. 
In this case, only the first reflexion contribute substan- 
tially to the image expansion (|16|l as if the second plated 
(at X = L) was pushed to infinity. The only remaining 
length scales are Xs and luj^ and the average force can be 
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approximated by 



KXs 



(24) 



with tA(M) = 3 /q cxp{-^/i 



cxp^— ^(7' + iM^j . This dimension- 
less screening factor is plotted in figure|Hl If the dynamic 
length Zjj^ is larger than the distance separating the shak- 
ing source from the plate, Xs/lui, 3> 1, the amplitude of 
the elastic field excitations are exponentially reduced be- 
fore reaching the plate's surface. More precisely, in this 
regime : 



exp 



i-V2u). 



(25) 



Oppositely. T a{xs/Iuj^) has a non vanishing finite 
limit if Xs/luj, goes to zero. then corresponds to 

the electrostatic force acting on a conducting plate 
facing a single "effective" charge 77s a'* located at 
(x = Xg, r = 0), so that, as in the quasi-static case above, 
F> (X -(7fea'^)V(Kxf-i). 

Fluctuations of the shaking-induced force: Beyond the 
definition and the characterization of the average inter- 
plate interaction, it is crucial to notice that the present 
shaking induced force is a time varying quantity. Rewrit- 
ing (|^ as F{t) = F^ [1 + cos(2a;st)], it is obvious that 
the force experienced by the immersed plates can not be 
a priori reduced to its sole average value. As a matter of 
fact the fluctuations of the instantaneous force are com- 
parable to its average amplitude. 



B. Homogenous Monochromatic Shaking. 

We now consider that the soft medium is uniformally 
shaken by rj^R, t) = 77s cos(cL'si). 77 is the superposition of 
in phase shakers of the -A- type uniformally distributed 
in space. All the elementary shakers produce periodic 
elastic deformations that "interfere" and give rise to the 
instantaneous forces F^{t) (resp. F^{t)) at x = 0"^ 
(resp. X = 0~) obtained from Eqs. (|2I6|) : 



d-l 



/ a,7^>^o(0,a;')dx' x 
Jo 

[l + cos{2ujst + 0>)] , 



(26) 



with 9' 



2arg 



!od.n>^^^oiO,x')dx' 



and 



F<{t) = -L 



d-l 



a,7^>,o(o,x')dT' 



[l + coa{2uJst + e<)] 



(27) 



with 6'^ = 2 arg [J^ d^TZ'^^ q(0, x')dx''\ . At every instant 
these two forces compete. They push the plate in op- 
posite directions to give rise to the net force F{t). The 



latter is exactly computed using E 
and after some elementary algebra 



F{t) = Fb 



The average Fb is given by: 



cos(2ijJst) 



[L/iV2LJ] 



(28) 



2 j-d-l 



Fb 



TBiu) 



viF 



B 



L 



1 cos(u/V2) 



2 cos(w/V2) + cosh(u/V2) ' 



(29) 
(30) 



The function Tb is plotted in section III-B, see Fig. 0] 
where its fast and slow asymptotics are qualitatively 
studied. 

Mean force tunability: The tunability of the present 
shaking-induced force is clearly revealed by Eqs. H29|l and 
%\[)\ . Precisely for fixed L: (i) The shaking pulsation lOs 
sets the sign of Fb to sign(cos[L/(v^^(.ijJ]), and its range 
to ~ l^^ (via the cosh term in (29)), (ii) then 77s allows to 
adjust the mean force amplitude arbitrarily. This second 
assertion is actually true only if i ^(^'^'^ with 
n integer. In such cases the plates experience no mean 
net force. 

Fluctuations of the shaking-induced force: Eq. H28|) 
indicates actually that F{t) "fluctuates" periodically 
around Fb with a pulsation 2ujs. Moreover the amplitude 
of the 2uOs components of the force is 1/ cos[L/ [l^^^/i)] 
times larger that its mean value Fb whatever the shaking 
parameters ujs and 77s. 

Shaking-induced organization: The two preceeding 
comments strongly suggest that the homogenous shaking 
may be used to dynamically induce spacial "localization" 
of the plates. As a matter of fact, when mesoscopic ob- 
jects are immersed in a soft medium, their intrinsic dy- 
namics (ignored in the present work) necessarily filters 
out high frequency effects of external forcing. If their 
response time is much larger than l/(2a;s), these "slow" 
moving objects only respond to the stationary compo- 
nent Fb of the shaking-induced force. The correspond- 
ing steady state interplate distance in the present simple 
geometry is then controllable by tuning Wg. 



C. Colored Noisy Shaking. 

The medium surrounding the plates is now supposed to 
be driven by stochastic shaking sources with zero mean 
and mean squared fluctuations given by Eq. H13|l . The 
noise averaged net force Fq experienced by the plate at 
x = is the sum of the mean values of the two indepen- 
dent forces F> and F<. Using Eqs. (12,18,20) this force 
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is written: 



Fc 



2L 



d-l 





[L 







{F>{x,Lu)+F<{-x,u))dx (31) 



Exact results: For a Lorentzian power spectrum, see 
Eq. H14() . the integration over frequency and position of 
all the elementary contributions Ffi^{x,uj) in l|31|l yields: 



Fc 
Tc (u) 



„2 drd-l 

rj^a i^ii 



^|^g(coth,-l) 

4~2 / coth (^V?^ + 



(32) 
(33) 



The scaling function Tc is plotted in Fig. [3 To achieve 
the exact computation of Fq, the integral over the in- 
plates q modes in (|33ll is performed. After some alge- 
bra we obtain: If c? = 1, Tc = (1 — u[coth(ii) — l])/u^, 
if d = 3, Tc = [3C(3) - tt^ - 6{u^ - ulog(l - e^") + 
3Li2(e^")]/(487rM^), where C is the Rieman Zcta function 
and Li2 the dilogarithm function : Li2(x) = X^fcLi /k"^ ■ 
Finally for the d = 2 case Tc = 2^JyK,^{2nu) + 



K2{2nu)], where Kn is the n K-Bessel function 

Notice that this last calculation can actually be by- 
passed, if one already knows the thermal Casimir force 
the left plate would feel if the soft medium had a finite 
correlation length ^, precisely if the elastic Hamiltonian 
were: H = («;/2)/[V(/)(i?)]2 + [4){R)/^fd'^R. Denoting 

this force Fcasimir (^^V^) i wc Can show that Eqs. H32l33f) 
can be recast into 1271: 



Fc 



2jkBT 



Fc, 



(0) - ^^Ca 



, (34) 



where i^casimir (0), of course corresponds to the Casimir 
force the plate experiences in a scale-free fluctuating 
medium, c.f. Eq. (|SJ). 

Qualitative approach: Eq. H31(l shows clearly that Fc is 
the result of the incoherent sum over frequencies and po- 
sitions of the forces produced by localized shakers of the 
-A- type. Note that the sum over the shaking frequen- 
cies is weighted by the noise power spectrum 5;^. This 
decomposition is now used to justify first the sign of the 
mean force, then the algebraic decays of the mean force 
as observed in Ea. H32|l and Fig. 

Only the elementary monochromatic shakers localized 
at distances \x\ < L on the left hand side of the plate 
have a challenger in the < a: < i inner region. Re- 
ferring to the discussion in Section IV-A, we know that 
> and that |i^^(a;,a;)| < |F^(-a;,w)|. Moreover, 
is necessarily a positive quantity for all frequencies. 
It then turns out that the net mean force is positive, 
for any choice of Suj: the two objects attract each other. 



Besides, F^[x,uj) ~ — F^(— x,^) for shaker located at 
x < L/2. This important property leads to the cancel- 
lation of the short distance divergences in the first term 
of equation (EH). Finally, \F^(x,bj)\ < \F^{-x,uj)\ if 
I a; I > L/2, so that Eq. H31|l can be approximated by 
Fc ^ 2Ll-\-^ I 1^1^ F< {~x, u:) dx. 

To sum up, the elastic field deformations produced by 
the ~A" type shakers located at x < —L/2 dominate the 
mean net force Fc. 

Without refering to any specific choice for S^, , the noise 
power spectrum is here characterized only by its width 
ils, see e.g. Eq. (13). Thus, going to the (x, l^j) variables 
and defining L* = max(F, ^fjj we easily deduce: 
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Fc 



^dL 



dx 



(35) 



since F^(x, lu) is exponentially weak if a; > and scales 
as l/x'^-i if X < c.f. Fig. (9). 

We can now integrate (|35|l and distinguish the two 
cases: 

(i) L/Iq^ 2> 1. In this wide noise spectrum limit Fc ~ 
iri^a'^n-'^/-f)L'l-^/L'^, the thermal Cas imir force scaling 
is obtained as anticipated in section III. 

(ii) L/Iq^ ^ 1. In this case the noise correlation time 
~ is larger than the field relaxation over distances 
smaller than L. The number of elementary shakers con- 
tributing to the force is substantially increased compared 
to (i). Depending on the space dimension d, we obtained 
the three scaling forms: if d = 3, Fc ~ {r/^a^ k~^)L'^^/ L , if 
d = 2, Fc - {Tia^K^^)L\\ log{lnjL), an finally if d = 1, 

^ (^78 oK~^)Zf2a • In agreement with the above exact 
calculations and with Fig. |31 

Fluctuations of the shaking-induced force: To conclude 
this section, we address briefly the question of the fluc- 
tuations of the force. In order to asses their relative 
importance with respect to the mean force value, we 
introduce the dimensionless ratio: Fc/AF, where the 
mean squared deviation AF^ is defined by: AF-^ = 
(F(i)^) — F^. Whereas Fc is the result of the compe- 
tition between the mean forces F^ and F^ of opposite 
signs, the forces' fluctuations on the two sides add up 
and AF^ = (AF>)2 -|- (AF<)2 with obvious notations. 
Note that the noise's four points correlation function is a 
priori required to compute AF. We make here one more 
assumption and consider that the r](R, t) are Gaussian 
random variables with correlations given by Eq. H13I14|I . 
This allows to fully characterize the shaking process and 
consequently the variance of the shaking-induced force. 
For sake of clarity we present without additional details 
the scaling form of the mean over variance ratio (ignoring 
logarithmic corrections): 

ik.('i^'^/('i:)Il ,36, 



AF 



L 
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we have also assumed d > 1, since the ID case deserve a 
more careful analysis as in the equilibrium context more 
carefully studied in |23|. The above expressions shows 
that, the present shaking induced force can not be a pri- 
ori reduced to its sole Fq value Fq- Indeed, depending on 
the geometrical aspect ratio L\\/L and on the noise cor- 
relation time rijT^ the variance of the force can dominate 
its mean value by orders of magnitude. 

D. Toward more complex dynamics 

To what extent do the results presented above apply to 
systems evolving according to more complex dynamics? 
It is noticeable that only the explicit form of the scaling 
functions Tx requires a precise description of the field dy- 
namics. Conversely, all the qualitative analysis, as well 
as the asymptotic expressions for the various forces, in- 
volve only the definition of the dynamical length l^^j. We 
thus expect that they should be extendable to any other 
dynamic scheme relating algebraically the spatial and the 
temporal relaxation scale. 

V. SUMMARY AND OUTLOOK 

To summarize, a simple model has allowed us to show 
that the concept of fluctuation-induced interactions can 
be extended to complex fluids driven to out of equilib- 
rium. Our main result concern identical plate-like in- 
clusions immersed in an homogenous medium externally 
driven by :-A- a monochromatic localized shaking source 
located between the plates, -B- a monochromatic and 
homogenous shaking and -C- a noisy colored shaking. 
These external processes generate effective forces on the 
plates without acting directly on them. The main fea- 
tures of the mean forces induced by these three simple 
shaking arc summarized in the table below and compared 
to their equilibrium thermal Casimir analog. 



Shaking 


Amplitude 


Range 


Sign 


A 






repulsion 


B 


ril 




tunable 


C 


v! 


power law 


attraction 


Thermal 


kBT 


power law 


attraction 



Five points are worth a highlight: 

(i) Whereas a shaking at a single frequency leads nec- 
essarily to short range interactions, noisy shakings 



generate effective forces that decays algebraically 
with the interplatc distance. 

(ii) If the external shaking occurs only between the two 

plates the resulting interaction is always repulsive. 
Conversely, when distributed on the whole sample 
the external shakers generate forces on both sides 
of the plates. The resulting net force is attractive 
for spatially uncorrelated shaking. On the contrary, 
with coherent shaking, interference phenomena be- 
tween the elastic medium distortions can lead to 
attractive or repulsive interactions. 

(iii) The soft medium-plates coupling has been here 
modeled by Dirichlet boundary conditions. As in 
the equilibrium case our results are expected to 
hold for other static boundary conditions represent- 
ing strong interactions. In contrast the description 
of situations of weaker coupling and/or cases where 
dynamic boundary conditions are required could 
yield to new phenomenology. 

(iv) The generalization of our results to an elastic 
medium characterized by an intrinsic relaxation 
scale ^ (see the discussion above Eq. (34)) can 
be done without much effort following blindly the 
procedures described in this paper. The only qual- 
itative difference in the results summarized in the 
above table is the forces' range: For shakings of the 
-A- and -B- type the range of the force becomes: 
min(^, l^). In the noisy shaking case -C- the force 
decay is not algebraic anymore but exponential for 
interplate distances larger than ^. 

(v) Beyond the derivation of their average values, we 

have also shown that the shaking-induced forces 
are in general expected to be strongly fluctuating 
quantities, whatever their precise origin. 

To conclude, let us recall that our description is clearly 
extremely simplified and of course is not supposed to 
model any real system. Even within this simple frame- 
work the extension of our result to many moving objects 
remains complex. But, the analysis of simplified situa- 
tions such as the one described here could certainly help 
understanding ordering and dynamic behavior of inclu- 
sions in out-of equilibrium complex fiuids. 
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